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Spectral sum rules for the quark-gluon plasma 

P. Romatschke and D. T. Son 
Institute for Nuclear Theory, University of Washington, 
Seattle, Washington 98195-1550, USA 

We derive sum rules involving the spectral density of the stress-energy tensor in 
= 4 super- Yang-Mills theory and pure Yang-Mills theory. The sum rules come 
I from the hydrodynamic behavior at small momenta and the conformal (in the case 

\ of = 4 SYM theory) or asymptotically free (as for the pure Yang-Mills theory) 

\ behavior at large momenta. These sum rules may help constrain QGP transport 

, ^ , coefficients obtained from lattice QCD. 

(M : I. INTRODUCTION 



port properties of the strongly 

m 

. Attempts have been made 
9|. These calculations rely on the 



i-^l Recently, much interest has been concentrated on the transj: 

coupled quark-gluon plasma (QGP) created at RHIC 1, 
to extract these coefficients from the lattice B, H, B B 

reconstruction of the real-time spectral function from Euclidean (imaginary-time) correlation 
functions, which for numerical data is an ill-defined procedure unless extra assumptions are 
^ \ made. In practice, the reconstruction amounts to postulating a form of the spectral density, 
\ and then fitting the parameters of the ansatz using lattice data. 

a^ 
rn 

cn 



Clearly, it would be of great help if some constraints on the spectral density can be 
derived. For nonrelativistic fiuids, there exist sum rules (for example, the /-sum rule) that 



O \ constrain the spectral densities 10| • One may wonder if such sum rules exist in a relativistic 
Q ! theory. Some progress has been made in this direction: for instance, Kharzeev and Tuchin 
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(and latter Karsch, Kharzeev, and Tuchin 12]) wrote down a sum rule for the spectral 
^ , density of the trace of the stress-energy tensor. The slope of this spectral density at zero 
frequency is the bulk viscosity Q. Although the sum rule does not fix the form of the spectral 
density, with some assumptions the authors of Refs. [ll|, [l2| argued that the bulk viscosity 
becomes large near the QCD phase transition. (As we shall see below, the precise form of 
our sum rule in the bulk channel is slightly different from that of Kharzeev and Tuchin, 
but some features of the later remain intact. We point out that the difference stems from 
a subtle non-commutativity of limits.) Several sum rules are also argued to hold for weakly 
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coupled relativistic theories by Teaney 

In this paper, we derive certain sum rules for the spectral density in hot gauge theories. 
We start with the A/" = 4 supersymmetric Yang-Mills theory, which is a prototype of the 
strongly coupled QGP, and derive the following spectral sum rule: 

2 _ 2 fduo 
5 -n J uj 

where p = — ImG/j(a;), and Gr is the retarded propagator of the T^^ component of the 



[P(^) -PT=oH], (1) 
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stress-energy tensor. Here pt=o{u!) is the spectral density at zero temperature, and e is the 
finite-temperature energy density. Besides AA = 4 SYM theory at infinite 't Hooft couphng, 
Eq. ([1]) is vahd for any theory whose gravitational dual has purely Einstein gravity. It is 
also valid for A/" = 4 SYM theory at any nonzero coupling. 

Another sum rule relates a linear combination of second-order hydrodynamic coefficients 
with the spectral function, 

oo 

1 2 [ duo 

2 Ti J UJ'^ 



where r] is the shear viscosity, and (the relaxation time) and k are defined in Ref. 15 
We then move to the bulk sector of QCD and show the following sum rule 

3(e + P)(l - 3c^) - 4(e - 3P) = - /- [p^^^'^u) - p^^olc^)], (3) 

where Cg is the speed of sound and p^^^^iu) is the spectral density for the trace of the energy- 
momentum tensor T^. Equation ([3]) is similar, but different, from the sum rule suggested 
by Kharzeev and Tuchin, and by Karsch, Kharzeev, and Tuchin. We shall show that the 
sum rule is a consequence of the hydrodynamic behavior of the QGP at large distances, and 
of asymptotic freedom at small distances. 

The structure of our paper is as follows. In Sec. [TTl we remind the reader how spectral 
sum rules can be derived. In Sec. Illll we derive two sum rules for the spectral function in the 
shear channel. One sum rule relates the total energy density with the spectral density, and 
another sum rule relate a linear combination of second-order hydrodynamic coefficients with 
the same spectral density. We verify both sum rules by numerically computing the spectral 
integral. We also comment on the possible form of the shear sum rule for pure Yang-Mills 
theory. In Sec. IIVI we turn to pure Yang- Mills theory and derive a sum rule for the bulk 
channel. 



II. KRAMERS-KRONIG RELATION 

For definiteness, consider the retarded correlator of the T^^ component of the energy- 
momentum tensor in, say, the A/" = 4 SYM theory. The spectral function is defined to 
coincide, up to a sign, to the imaginary part of the retarded Green function of T^^, 

Pi^, q) = - Im Gr{uj, q), (4) 

where we assume q to be along the z direction. Since Gr is the Fourier transform of a real 
function (recall that Or determines a linear response), we have 



Gif(-u;,q) = G*ji{uJ, -q). 



(5) 



FIG. 1: Integration contour for the integral 



Let us consider the function fq{Z), defined so that /q(c<J^) = 6*^(0;, q), which has a cut 
from Z = to Z = 00. 

We can write down a Kramers- Kronig relation for /q(^). Pretend for a moment that 
fq{Z) — > as Z — * CX3. Taking the integral of fq{Z)/ {Z + a^) over the contour in Fig. [1], we 
find 

if fq{Z) does not have singularities except for the positive real semi-axis. In a relativistic 
field theory, typically Gr diverges as ^ 00. For example, the T^^ correlator Gr{u) grows 
like Lo'^lnu at large u (see, e.g., One can subtract this zero-temperature piece G'^"°(u;) 

and we denote SGr = Gr — G'^^. But, as will be shown, there remains a constant piece 
which also needs to be subtracted. Therefore, we define 

Uuj^)=5Gr{uj,c^)-5G^, 5G^= hm 5G^(cu,q), (7) 

OJ—t+tOO 

so that the Kramers-Kronig relation ([6]) for this function fq{Z) is valid. 

In a conformal, large- iV theory^, the asymptotics of G'_r(c<j) at small u and k is known to 
second order from hydrodynamics 
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Gr{uj, k) = Gr{0) - irju + {r]T^ - u;^ - ^nq^ + .... (8) 

The constant Gr{0) may depend on the way the correlator is defined. When the correlator 
is defined through the response to metric perturbations, Gr{0) = P. For this particular 
correlator, fq{Z) does not have singularities outside the positive real semi-axis and the 



^ Outside the large- limit, there can be non-analytic terms (such as w'^/^) present in Gr{u)) (141] . 
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contour integral may be deformed to give 

oo 

f g{-a') = SG r{0) - 6G^ + r]a+{r]T^-^) a' -'^f^q' + ... = -- [ du p{u , q) , (9) 



where 5p{uj, q) = p{uj, q) - pT=o{'^, q)- 
Setting a = in this formula, we find 



oo 

- /,(0) = 5G^ - 6GniO) + \Kq' + 0{q') = - 6p{uj, q), (10) 



which for q = will become Eq. ([T]). Subtracting the a-independent part, the Kramers- 
Kronig relation becomes 



oo 



12, 2a^ fdu 6p{u,q) 
T]a + {7]T^ - ^Kja H = / 2 I 2 • (11) 



TT J UJ 00"^ + a 

From this we derive another sum rule: 



oo 



r/r. - 1/. + 0{q') = - f—, [5p{uj, q) - ^ , (12) 

71 J UJ^ 







which for g = is Eq. ([2]). Note that the coefficient k enters this equation, which is 



15 



interesting because in the hydrodynamic equations k couples only to curvature tensors 
and hence drops out for flat space. Equation (fT2!) suggests that /t can be determined from 
fiat-space physics. In fact, k, can already be obtained from Euclidean correlators at = 
and small q. For weakly coupled SU(A^) gauge theory, we find a nonzero value for n at the 
lowest order of perturbation theory (see Appendix Curiously, n divided by the entropy 
density only differs by a factor of about two between strongly coupled A/" = 4 SYM and free 
SU(A^) gauge theory. 



A note on the definition of the correlators 

It is clear from the previous discussion that, in order to derive the sum rule in a particular 
theory, one should use the same definition for the correlation function in the UV {uj oo) 
and IR (cj — > 0). In this paper, we use define the correlators through the partition function 
Z in curve spacetime. The one- and two-point functions are given as the first and second 
derivatives of \nZ with respect to the metric. In Euclidean signature, one has 

5\nZ = l [dx {T>^''ix))5g^,{x) + l [ dx dy {T^^''ix)Tf"^{y))5g^,{x)5g,M + ■ ■ ■ (13) 
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In other words, 



52 InZ 



An alternative definition of tlie correlator is through the path integral in flat space, 

1 



{T^^{x)T'"^{y))' 



VAe-^^T^''{x)T'"'{y), 



(14) 



(15) 



where A represents all fields in the theory and Se is the Euclidean action. The two correlators 
differ by a contact term^ 



(T^'^(x)T^"(l/)) = {T^'^{x)TP''{y)y -A 



6^Sf 



(16) 



.Sgf,uix)6gp„{y) 

Analogously, one can define the Minkowski-space correlation functions. The retarded Green 
function is found from the linear response, 

5 



{T^^ix)TP"iy)), 



- 2 



(17) 



The advantage of using the correlator defined through Eqs. f|T3|) and f|T7|) is that we know 
this correlator at low momenta through hydrodynamics. Indeed, using the hydrodynamic 
equations one can establish how a system responds to external gravitational perturbations, 
and then use Eq. (HM to find the correlation functions (see e.g. 17|, ll8|]). In addition, this is 
the most natural definition that comes out of AdS/CFT correspondence. Note that on the 
lattice, so far what is normally measured is (fT5|) . However, this difference does not matter as 
far as the Kramers- Kronig relation (Q is concerned, because subtracting SC^ from 6Gfi{0) 
makes the contact term drop out from fq{0). 



III. SHEAR SUM RULES IN 7^ = 4 SYM AND PURE YANG-MILLS THEORY 



A. AdS/CFT Calculation of /,(0) 



For the case of large 't Hooft coupling, properties of A/" = 4 SYM can be calculated 



19 



20 



2l|. In particular, it is known how to calculate finite 
To find the {T^^T^y) correlator, we 



using the AdS/CFT duality 

temperature correlators Gj?(u;,q) in AdS/CFT |l7l 
solve the equation of motion for the xy component of the metric, which is essentially the 
equation for a minimally coupled scalar, 

1 



-9 



-99 



9^ kpkp 







(18) 



^ A contact term is a term in {T^^ {x)TP'^ [y)) that is proportional to 5'^{x — y) or its derivatives (corre- 
sponding to a constant or polynomial in momentum space). For theories which do not have derivatives 
of the metric in the action (such as Yang-Mills theory), the contact terms can only be constants. 
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where /i = ... 3 indexes the usual four field theory dimensions and g is the determinant 
of the (five-dimensional) metric. We denote the fifth dimension by z (not to be confused 
with the spatial direction in the previous section), where z = corresponds to the four 
dimensional boundary of AdSs space. Finite temperature correlators in AdS/CFT can be 
studied by considering the metric of a static black hole in the bulk. The location of the event 
horizon zh of the black hole is related to its Hawking temperature, zh = (ttT)"^. Being 
interested in SGr at large imaginary u; ^ T, we can restrict ourselves to the region of AdS 
space very close to the boundary. 

For convenience, we shall use here the metric in Fefferman- Graham coordinates, which 
has the following form [23] near the boundary at z = 0, 

_//y-2 _|_ _|_ J 2 fi'2 2 

^ ^2 at +(^^ ^ ^ (3^^2 ^ ^^2) ^ Q^^4) ^ (ig) 

z AZfj 

where R is the scale set by the AdS radius. In these coordinates, Eq. (!T8!) becomes 



+|l + ^l-V-(l-J^)qV = 0. (20) 



We will put q = from now on, and consider Euclidean momentum tu^ = —Q^. In the deep 
Euclidean region +oo, most of the interesting dynamics happens near the boundary, 

so can solve Eq. fl20|) iteratively in inverse powers of zh- We expand the solution as 

= 00 + 01 + ■ ■ ■ , (21) 

where 0o = \{QzYK2{Qz) is obtained by sending zh oo in Eq. ( !20l) and demanding 
regularity at z — > cxd (see also [l3|). The first correction 0i satisfies 

0'i' - -0'i - Q'0i = 3 = 3Q^-^0o . (22) 

Z ^Zjj 

The solution to this equation is formally given by a Green's function, 

0i(z)= jdz'G{z,z')j{z'), (23) 

where G{z,z') can be constructed from the known solutions to the homogeneous Eq. (122]) . 

h{z) = {QzfK^iQz) , h{z) = {QzfUQz) , (24) 

as 

= nf J fii^)Mm^ - + f2{z)h{z')e{z' - z)] , (25) 

where the Wronskian W[fi, /a] = /1/2 - f[f2 evaluates to W[fi, f2]{z') = Q^z'^. 
Evaluating the integral fl25|) using 

00 

ldzz'Kl{z) = ^, (26) 
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we then find the small z asymptotics of - 



1 to be given by 



^i(^) = -T7T7r- 



Recalling that the correlation function is given by 16 1 



which at T = reproduces the well known result 

iV2 4>',{^z) 



G 



T=0, 
R ' 



32^ 



rQ' In Q 



The first correction due to temperature is 



lim dGftiuo) 



3 

20^4" 



20 



2rr4 



(27) 



(28) 



(29) 



(30) 



There are also contributions to the correlators from the boundary terms in the action (i.e., 
J d'^x a/— 7, cf. ji^l) but one can check they contribute the same amount at any u, and they 
are the only contribution at = (the contact terms). Thus, in iV = 4 SYM theory one 
has /g=o(0) = — |e and hence the sum rule (ITOl) becomes 



duj 



UJ 



[p(^) - Pt=o(^)] • 



(31) 



We will call this sum rule the "shear sum rule" as the slope of p{u!) at = is the shear 
viscosity. It is clear from our derivation that this shear sum rule holds in any theory with 
an Einstein gravitational dual. 



B. Rederivation of /g(0) in Af = A SYM from OPE 

Within Af = 4 SYM theory the sum rule can be derived without relying on gravity. We 
start with the operator product expansion (OPE) of the stress energy tensor 24 



T,.{x)TpAO) ~ + i,.p.„/5(x)T„^(0) + ■ ■ • (32) 

Here A contains various Lorentz structures, all scaling as x~^, and is given explicitly in 



Ref. [2J] in terms of three constants a, b, and c. In a thermal ensemble the second term 
averages to a constant contribution to the correlator. Setting fi = p = x,h' = a = y, and 
performing Fourier transform, we get^ 

6Gii{uj) U^ioc = ]^(^±3_p . (33) 

Kv ;n *oo 14a -26 -5c ^ ^ 



In Ref. [24!\ the correlator is defined as the second derivative of the partition function with respect to the 
upper components of the metric, while we differentiate In Z with respect to the lower components. This 
difference, together with the sign change by going from Euclidean to retarded propagator, has been taken 
into account in Eq. 
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In = 4 SYM theory the coefficients a, b, and c are given by 



25 



16 



9n^ 



{Nl - 1), 



17 



97r6 



{Nl - 1), 



92 



97r6 



and hence 

On the other hand, the shift of Gr when uo 
predicts 5Gr{Q) = 



11 



P. 



(34) 



(35) 
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can be found from hydrodynamics, which 
17|. Therefore we find /g=o(0) = — f-P; which corresponds to ([T]) 



since for a conformal field theory P = ^e. 

A remark is in order. It is known that the constants a, b, and c are independent of the 
coupling in A/" = 4 SYM theory: in fact, their value can be found from one- loop calcula- 
tions 



25 



Therefore, the sum rule is valid for any nonzero value of the coupling. 



C. Calculation of /g(0) in pure Yang-Mills theory 

In pure Yang-Mills theory, the UV behavior is that of a weakly coupled field theory. The 
leading terms in the OPE are the same as for free fields. The coefficients a, b, c can be found 
from the general formulas 24 

^ n^- -T^^v, (36a) 



277r6 ^ TT' 

4 1 



118 , , 

-n^ ^nf -n^, (36c) 



where Ug, nj, and are the number of real scalars, Dirac fermions, and gauge fields in the 
theory. 

For pure Yang-Mills theory, by repeating the calculations in section IIIIBt we find 
fq=o{0) = —2P (this can be checked directly by computing the relevant Feynman diagram; 
see Appendix [B]). However, there is an additional subtlety here that was not present in 
the previous section because the OPE of two components of the stress-energy tensor may 
involve terms like a^F^, where as is the strong coupling constant. Though formally higher 
order in a^, these terms average to e — 3P which is a constant independent of the scale x in 
Eq. fl32|) . Therefore, we can tentatively write a sum rule for pure Yang-Mills theory. 



CO 



'-^ + C{e-3P) = - [— [p{u) - pT=o(^)] . (37) 



where the constant C is left to be determined in by a more accurate calculation. (How- 



ever, naively applying the results of Refs. [26| about the absence of the leading order gluon 



condensate contribution in the tensor glueball channel would imply that C = 0). 

In the large- A'c limit, where second-order hydrodynamic coefficients are well defined, the 
sum rule ffT2l) is valid (except that one has to use a proper definition of k for nonconformal 
theories [27|). 
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FIG. 2: Numerical results for the spectral function Sp{uj, 0) for strongly coupled AA = 4 SYM. The 
results shown correspond to the integrand of the sum rules ([I]),©, respectively. The inset in (a) 
demonstrates the near-exponential drop in the amplitude of Sp{uj,0). Horizontal lines are visual 
aids to the eye. 



D. Numerical verification of sum rules in AdS/CFT 

In AdS/CFT, the spectral function p{uj) can be calculated numerically for arbitrary 
frequency /moment a from the solution to the mode equation ( fT8l) . For convenience, we 
adopt a metric and coordinates such that Eq. ( JTSl) becomes 16 



uf{u) 



>>' + 



w 



up{u) 



0, 



(38) 



where u = z^/zj^, f{u) = 1 — m^, w = Lj/{2nT), and q = g/(27rT). Again, we set the spatial 
momentum q = in the following, although the method described below can also handle 
non- vanishing momenta. 



We follow the algorithm by Teaney 13|], which is outlined here for completeness, fixing 
some typos in Ref. [l3|^. Equation (!38|) is recast in a system of coupled first-order equations 
by introducing vr = 0'. Discretizing derivatives as 0' = + 6u) — (p{u)) /6u, Eq. (!38|) can 
be explicitly integrated forward from a point close to the boundary u = Uq. If vr is taken to 
be defined at half-integer stepsizes t:{u) = 7r((n + l/2)6u) and (p at integer stepsizes (f){u) = 
7r{nSu), then the resulting algorithm is second-order accurate in 6u ("leapfrog algorithm"), 
suggesting numerical stability. To start the algorithm, initial conditions for and vr at 
u = Uq, u = Uq + ^Su, respectively, need to be specified. For uq sufficiently close to the 



A version of the C++ code will be made available at http://hep.itp.tuwien.ac.at/~paulroin 
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boundary u = 0, Eq. (138|) may be solved analytically, yielding the pair of solutions 

,2 



W 



2 

$2(m) = \og{u)^i{u) + 1 + W^M w^u^ + . . . , (39) 

2 9 

where the u"^ term in $2(w) is arbitrary (it can be any multiple of $i(m)) and was set to 
zero in accordance with the convention by Kovtun and Starinets 28|]. The analytic result 
for ^i{u) (and its derivative) is used as initial conditions for 0(Mo),7r(Mo + ^Su), which can 
then be integrated forward to give a numerical solution 0i(ui) with ui close to the horizon 
u = 1 (the same procedure for $2('w) gives 02 (""i))- The physically interesting solution for 
0(ti) is the one that corresponds to an incoming wave at the horizon, (j){u) ~ (1 — m)~™/^. 
Solving Eq. fl38|) analytically close to u = 1 one finds for the incoming wave solution 



,2«w^ + 3w^-?w ,,w(4w^ + 7?.w2 - 2w + 4z) 

1 - (1-u) -, (1-ur— — ^ + 

^ ' 4 1 + w2 ^ ' 32w2 + 3iw-2 



(40) 

The real solutions 0i(m) and 02 (m) are linear combinations of the incoming and outgoing 
wave solutions, 

0i(w) = A(w) r\u) + i?(w) r'{u) , 

02(m) = C(w) 0^°^(n) + D(w) 0''^=(m), 

where 0™'^ denotes the complex conjugate of 0™'^. For given w, the complex constants 
A, B, C, D are, e.g., calculated from the numerical solution of 0i^2 at u = ui and u = Ui — 5u 
and the analytic solution (HOl) for 0™"^ close to the boundary. As a consequence, one can 
construct a numerical solution to Eq. ( l38l) with incoming wave boundary conditions by 



^incnum/ N _ ^ ^ / N ^ ^ / N /42) 

^ ^ ^ AD-BC^^ ' AD-BC^^^ ^ ^ 

which can be normalized to 0(0) = 1 by realizing $i(0) = 0, $2(0) = 1, so that close to the 
boundary 

0i-"(uo) = -^'^•i(«) + . (43) 

In practice, we found the choices Uq = 10~^, Ui = 0.999 to give acceptable numerical accu- 
racy. Once the normalized solution to the mode equation is known, the retarded correlator 
is obtained from 

TT^N^T^ a,0'-™(no) 
Giiiuj, 0) = hm ^ . (44) 

4 uo-*0 Uq 

In particular, using the analytical results for $12 one finds for the spectral function 

Sp{u;,0) = -4P^-2nPw\ (45) 
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where P = ^N'^T^ in strongly coupled = 4 SYM. The numerical result for the spectral 
function is shown in Fig. [21 As can be seen from this figure, 5p{u!, 0) first increases as 
a function of u, reaching a maximum at around w = 0.45, then decreases strongly and 
oscillates around zero with an amplitude that decays quasi-exponentially. Fig. EJ^b) shows 
the spectral function where the leading hydrodynamic behavior rju has been subtracted. As 
can be seen from this figure, for small frequencies the spectral function seems to behave as 



r]uj 



aow^ + aiw^ + O(w^) , (46) 



where numerically we determine ao — 1.72, ai ~ —3.0. 

With the value g^/q=o(0) = ^ calculated above, the first sum rule ([1]) would imply the 
identity 



i=o.i5=ir^f-^-v 

20 TT./o w V B(w) 2 



1 Z"'^'""'' dw ( D(w) TT , , 

-w"^ = 0.1500008(44), 



(47) 



TT Jq w \ B(w) 2 
which we can confirm up to five digit accuracy when choosing Wmax = 6 in practice. For the 

1^ _ op l-log(2 
2'^ (27rT)2 



second sum rule ([2]), r^r^r — \k = 2P ^^'"rla^ from Ref. [15| imphes 



1 - log(2) . 0.306853 I T '-^ ("^ " "^^ " "1 

' VT Jo w3 V 5(w) 2 2 J 



, dw ( Diw) TT 4 W ■ 
+ - / — -^-^ - TTw' - 77 = 0.30686(2), 



TTWmax TT J W"^ \ B{w) 2 2 



(48) 



indicating that the numerical result matches with four digit accuracy. While it is possible 
to improve the numerical accuracy further, we take this agreement of at least one part in 
10~^ between the analytical and numerical results as an indication that for A/" = 4 SYM, 
the sum rules (fTl) and (12]) are correct. 



IV. THE BULK SUM RULE IN QCD 

In this section we revisit the sum rule satisfied by the imaginary part of the correlation 
function of the trace of the stress-energy tensor T^. As the spectral density in this channel 
is related to the bulk viscosity, this sum rule will be called the "bulk sum rule." We show 
that this sum rule indeed exists, but its form is slightly different from the one given in 
Refs. HQ. 

In this section we shall be concerned with metric perturbations of the following form 



(49) 



12 



or (5o 



-20, 



1), with f2 ^ 1. For these perturbations, the partition function 



expansion defined in Eq. (fT3!) can be exphcitly given as 



6lnZ 



+ - I dxdyr^^,{x)Tip„{y){T^''{x)TP''{y))n{x)n{y) . (50) 



We then define the correlators of 0{x) as follows: 
5\nZ 



{0{x)) 

{e{x)e{y)) 



5VL{x 
6^ InZ 



Q=0 



6n{x)6n{y) 
6 



sn{x] 



~ggAT"'{mn=o = V,.{x){T^''{x)) , (51) 

(52) 



where we recall that in conformal field theories {0{x)) = {6{x)6{y)) = 0. As a consequence of 
the definition (1511) . the correlator {6{x)6{y)) differs from Tjniyfipfj {Tf^'TP") by a contact term, 
which in its turn differs from (Tj^T^)' by a contact term. Our subsequent calculations are 
simplest when using the correlator {6{x)6{y)) defined in this fashion. 

Consider the pure Yang-Mills theory. We need to know how to couple Yang-Mills to 
an external metric perturbation of the form (H9|) . This is done through changing the bare 
coupling g'^ = iiras, so that it is dependent upon the metric. In particular, the Euclidean 
action of pure Yang-Mills becomes 



dx 



where we have rescaled the gauge fields so that the field strength tensor is given by 



and / are the SU (N) structure constants. From this we find 



6Sf 



Pigs 



dS 



Pigs) p2 



dgs 

where Pigs) = ^d^gs is the beta function. As a consequence, one has 



(53) 



(54) 



(55) 



(Oix)) 



Pigs 
2g's 



{FUx)) , 



(56) 



{dix)d{y)) 



Pigs 
2g's 



{F\x)F\y)) + (3igs) 



d (Pigs) 



dgs V 2^, 



{F')5ix - y) . (57) 



13 



If we are interested in computing Gr{(jj), it is most convenient to choose A ~ cu, so the two 
terms in Eq. (157|) can be evaluated perturbatively without large logarithms. In the weak 
coupling regime 

(^{gs) = -hogl-hgl + --- , (58) 

so the first term in Eq. (!57ll is proportional to gj{uj)T'^, while the second term is proportional 
to g^{uj){e — 3P). When u oo, the correlation function vanishes because of asymptotic 
freedom. 

Now let us compute (06)^ at small frequencies. For that we need to find the response of 
the system of an external metric perturbation with Q = Q{t) varying slowly with t. Since 
the perturbation is spatially homogeneous, we expect the fluid to remain at rest (m° = e^, 
ti* = 0), but the temperature will have time dependence: T = T{t). It is more convenient to 
work with the entropy density s instead of T. When metric perturbations are slow, entropy 
is conserved. The solution to the equation for entropy conservation, V^{su^) = 0, is 

s = e^^so , (59) 

and therefore 

This result may now be directly applied to the definition of the correlation function, and 
hence we find 

{99Uu - 0, q = 0) = ^{V^rp = - (3.^ - 4) (e - 3P) . (61) 
Now let us derive the spectral sum rule. Introducing the spectral function 

^buik(^^) in the 

bulk channel, we find 

00 

3,|_4)(.-3P)^1/^*P-M. (62) 



where e and P are now the thermal parts of the energy and pressure (with the divergent 
vacuum contributions subtracted out). This, we argue, is the correct version of the sum 
rule by Karsch, Kharzeev and Tuchin [ll|, [l3]- The right hand side can be transformed into 
Eq. ([3]) by using the thermodynamic relations de = Tds and dP/de = c^. Note that Eq. (!62|) 
does not coincide with Refs. {ill. E^. which had Td/dT instead of 3sd/ds. The issue is the 
non-commutativity of the g — and u —* limits in the bulk channel (see Appendix O) . 
The correct expression for the right hand side follows directly from entropy conservation in 
hydrodynamics. 

In the weak coupling limit of high-temperature gauge theory, the pressure is given by 



30 



P = T'{A + glB + 0{gl)) , (63) 
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where A, B are constant that are unimportant for the following discussion. Calculating the 
trace anomaly from e = — P one finds 



e - 3P = 2BT%(3{g,) ~ 0{gtT 



which implies 



31 



3.- 

ds 



4 ) (e - 3P) ~ Oig^T') 



(64) 



(65) 



where ^ = 3 + P{gs) was used. On the other hand, the integral over the spectral 

In order for our sum rule (1^2]) 



function gives a contribution O^g'^T^) at low frequency [31 
to hold for weakly coupled QCD, the 0{g'^T'^) contribution must be canceled to leading- 
order by the high frequency tail in the spectral function. There are indications that this is 
indeed what is happening in weakly coupled QCD 



36 . 



V. CONCLUSION 

In this paper we have written down several sum rules involving the spectral functions in 
hot gauge theories. The sum rules can be checked in A/" = 4 SYM theory using gauge/gravity 
duality. The bulk sum rule for QCD was also derived. We still have some uncertainty in the 
shear sum rule in QCD, but hopefully this will be resolved in the future. 

Some conclusions may be drawn from our work. First, note that the left-hand side (LHS) 
of the sum rule is positive in strongly coupled Af = 4 SYM because ?7r^ > l'^- This 
implies that the spectral function 6p in the shear channel must be larger than r]uj for some 
frequencies, or otherwise the integral would not be positive. This feature of the spectral 
function is clearly seen on Fig. [2l This requirement is not satisfied by the simplest Lorentzian 
ansatz for the spectral function, Sp{uj) ~ ^jr^- A similar argument can be made for weakly 



coupled QCD in the large Nc limit, because there ~ Qri/{sT) 32] and rj/s ^ a^^lna^^^^ 
33|, while n T'^ (see appendix so one expects r^r^r > Therefore, to satisfy both 
shear sum rules in QCD, an ansatz more sophisticated than the simplest Lorentzian ansatz 
6p{uj) ~ ^2%2 is needed. 

Moreover, the LHS of our sum rule for the bulk sector can be evaluated using lattice 
results for the thermodynamics [35?]. The result turns out to be negative for all temperatures 
above the deconfinement transition (see Fig. [3]) Fig. [3] demonstrates that the sum rule ([3]), 
cannot be directly used to extract information about the value of the bulk viscosity in QCD, 
unless additional phenomenological assumptions are made, for example as in Refs. HQ- 
In our language, Refs. assume pt=o{^) to contain a "non-perturbative" part (asso- 

ciated with the phenomenological gluon condensate), that — once subtracted — is offsetting 
the negative LHS of the bulk sum rule Q. 



^ For non-conformal theories, corrections to the definition of k should be suppressed by an additional power 
of [34|. 
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FIG. 3: The LHS of the bulk sum rule [Eqs. ([3]) and (I62p ] (full line), as evaluated from SU(3) lattice 
data [35]. The result is negative for all temperatures shown, and should not be associated with the 
value of the bulk viscosity in SU(3). For comparison, the result from evaluating Eq. (jCTh (corre- 



sponding to the sum rule by Kharzeev and Tuchin 11[) is shown (dashed lines). The horizontal 
line is a visual aid to the eye. 

While this work is being completed, we became aware of related work by Simon Caron- 



Huot [36|. We thank L. Yaffe for discussions, G. Moore and S. Caron-Huot for comments on 
the manuscript, J. Engels for providing the lattice QCD data used in Fig. [31 and S. Caron- 
Huot for correcting an error in our earlier treatment of the bulk sum rule and for giving us 



a preview of Ref. [36|. This work is supported, in part, by U.S. DOE grant No. DE-FG02- 
00ER41132. 



APPENDIX A: THE COEFFICIENT k 

Here we calculate directly the Euclidean correlator 

G{xi,X2) = {T:^y{xi)T^y{x2))j, (Al) 
for a free SU(N) gauge theory at finite temperature T. We have 

G{x^,X2) = {{d.A: - d^Al) {d^Al - dyA:) - d,At) [d,Al - dyA^) )^ , (A2) 

where each of the building blocks is a correlator of the form 

C'iii2jij2hl2mim2{Xl} ^2) = {diiA^dj-^A-j^dlj^Af^dmiAfj ^^j, 

lyr ^i{Pi+P2)-xi+i{P-i+PA)-X2 ph ph ph pmi 

JPl,P2,P3,P4 

{^^{Pi)K'{P2)A[\P,)A^\P,))^ , (A3) 



X 
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and P = (p,P4) = (p, 27rTri), = TJ2nI (2^- ^^^^ interested in the Fourier 
transform of these correlators, 



«l*2ili2'l'2"lim2 



dr I d 



ili2jl32lll2m\m2 



[x, 0) 



X 



where 



P2,P3,P4 

+ {Al^Q - P2)A[^{P,))^ {Ai^{P2)Ar{P,))r 
+ {Al'iQ - P2)K%P2))t {A[\P,)Ar{P^))T] 

mPM'dP2))T = ;^5n,+n„o(2vr)=^5(pi + P2)KtiPi) , 



(A4) 



(A5) 



and A|^^(Pi) is the gluon propagator. Note that the last term in Eq. (lA4p corresponds to 
a disconnected diagram; we are only interested in the connected Green's function, so this 
term will be dropped in the following. Using Feynman gauge A^^(Pi) = dab^^" Pi^-, Eq- (1A4p 
simplifies to 

C^,^2n32hl2^.m2{ci,q,) = {N^ - l)^ {Q + K)-^ K'^Q + K^^ X 

[{Q + + K'^S'^^'^IQ + Ky^^^^^] , (A6) 

and hence the Fourier transform of the energy momentum correlator becomes 

G{Q) = {Nl - 1)^ (g + KY^K-^ X 

Jk (A7) 

[Aklkl -2K-{Q + K){kl + kl) + K^kl + (g + Kfkl + {K ■ {Q + K)f] . 

Since we are interested here in the case for vanishing external frequency {p^ = 0), the 
thermal sums are readily evaluated. Dropping the vacuum part and using the substitution 
k — > k — q in part of the integrand one finds 

T^(g + ir)-2 - ^^^^ 

n 

tY^{q + k)-^k-' 

n 

2^ in 



{K-{Q + K)f 



(g + KYK^ 




(A8) 



Expanding the integrand to (9(q^), all the remaining integrals can be done analytically and 
one finds 



5G(0, q) 
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(A9) 
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so that 



Ar2 



(AlO) 



APPENDIX B: THE COEFFICIENT /^(O) IN WEAKLY COUPLED SV{N) 

Calculation of = linig^^oo Gr{Q) starts similar to the calculation for k in the previous 
section, leading to Eq. ( ]A7I) for Gr{Q). For q = the sum-integrals become 

l + 2n{k) 



n 

n 

2^ (D 



2k 

l + 2n{k) 
2k 

l + 2n{k) 



1 



+ 



1 



g| + 2ikqi q\ — 2ikq4 
ikq4 ikq^ 



{K-{Q + K)f 
{Q + KfK^ 



2k 

l + 2n{k) 
2k 



g| + 2iA;g4 g| — 2ikqi 
Pqj k'^qj 



+ 



_g| + 2ikq4 g| - 2ikq4 



(Bl) 



where we used n{iq4) = n{2TTiTn) = 1. Evaluation of the remaining integrals is straightfor- 
ward and we find 

lim G(g4, 0) - G(g4, 0)t=o = fg=o{^) = -2P, (B2) 

(74^00 

where P = ^{N^ - l)7v^T^. This result can also be obtained by integrating the result for 
the spectral function from Ref. |37|. 

APPENDIX C: NON-COMMUTATIVITY OF THE u^O AND q^O LIMITS OF 

THE (66) CORRELATOR 



In Sec. IIVI we have shown that 

lim lim(^^)^(^,q) = 

Using the same method, we now show that 



hm \im {ee)j^{uj, q) 

This result is consistent with previous results 
tegral following the method of Ref. 



3^-4) (e-3P). 



(CI) 



(C2) 



38 



39 



derived through the Euclidean path in- 



It should be expected: Euclidean correlators are 
defined with discrete Matsubara frequencies ue = 2imT , and the only sensible zero momen- 
tum limit in the Matsubara formalism is to set uje = first, and then take q — > 0. 
We turn on a static metric perturbation. 



-2n{x) 



(C3) 
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When ri(x) varies smoothly, one can use hydrodynamics to find out the response. For static 



perturbations we expect the response will be static. The velocity field is m° 



u 



0, 



and the temperature depends on space, T = T(x). Substituting Tf^^ 
into the equation V^T'^'^ = 0, we find 

d,P - (e + P)din = 0. 

Using dP = sdT and e + P = Ts, the solution to this equation is 

T = Toe^. 

The correlator is found from 

{9e){uj = 0,q^ 0) 



(C4) 
(C5) 



d_ 



'3T11) = -{t^-^ (e-3P). 



(C6) 



One apparent paradox is that if one writes down the bulk sum rule for any spatial 
momentum q 7^ 0, the integral should be equal to — (6*6') (0,q) which is given by flC2p 
but not dCT]) : 



nn d 

T — 

dT 



4 



dui 



3P) = Z l-Sp''-'\u,q), q^O. (C7) 



There is no contradiction, however, as the integral in Eq. flC7p is expected to receive a finite 
contribution from the region ~ g, in particular from the sound-wave peak at uj = CgQ, as 
(96) correlator has a sound-wave pole. When g — 0, this region shrinks to zero size, but its 
contribution remains finite. The contribution from the sound-wave peak can be calculated 
as follows: to leading order in hydrodynamic fiuctuations, T^^ = T^^ = T^^ = CgT°°, so 



bulk 



R ■ 



(C8) 



Defining {T^'^T^^) r =^xduj . q) 
corresponding to Xl 
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CO 



^^^'^''^'^ and using Teaney's result for the spectral density 

r,gV2 



c.g)2 + (r,gV2f 



+ [00 



(C9) 



where Tg = {e + P) ^ [^r] + (^), the integral over the sound-wave pole at positive frequency 

00 

gives limq^o ^ Spl{uj, q) = e + P. Therefore, the contribution for the integral over 5p 



bulk 



is precisely the difference between the LHS of Kharzeev-Tuchin's and our sum rule. 



T 



d_ 
df 



3s 



d_ 

ds 



(e - 3P) 



3c' 



,2\2 



P)- 



(CIO) 



On the other hand, in the bulk sum rule fl62l) . we first compute the spectral function p{uj) at 
any finite, nonzero uj by setting q = in p(a;,q), and then take the spectral integral. The 
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sound-wave contribution does not appear in this integral, which means that our Eq. (l62l) . 
but not the Kharzeev and Tuchin's version, apphes. Note that at zero temperature the two 
hmits u ^ and q ^ commute. 
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